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Kalata Bl has been demonstrated to have bioactivity relating to membrane disruption. In this study, we 
conducted coarse-grained molecular dynamics simulations to gain further insight into kBl bioactivity. The 
simulations were performed at various concentrations of kBl to capture the overall progression of its 
activity. Two configurations of kBl oligomers, termed tower-like and wall-like clusters, were detected. The 
conjugation between the wall-like oligomers resulted in the formation of a ring-like hollow in the kBl cluster 
on the membrane surface. Our results indicated that the molecules of kBl were trapped at the 
membrane-water interface. The interfacial membrane binding of kBl induced a positive membrane 
curvature, and the lipids were eventually extracted from the membrane through the kBl ring-like hollow 
into the space inside the kBl cluster. These findings provide an alternative view of the mechanism of kBl 
bioactivity that corresponds with the concept of an interfacial bioactivity model. 

Cyclotides are 28-37 amino acids peptides found in various plants, such as Rubiaceae, Violaceae, 
Cucurbitaceae and Apocynaceae 1 . These peptides are classified into the Mobius, Bracelet and Trypsin 
inhibitor subfamilies 1 . All cyclotides contain a cyclic cystine knot (CCK) motif, a defining characteristic 
that is a combination of a head-to-tail cyclic backbone and a knot of three disulfide bonds that connect six 
conserved cysteine residues 24 . The amino acid sequences between each pair of the conserved cysteines are termed 
loop 1-6 2 . Cyclotides are stable in various physical and chemical conditions as a result of the CCK motif' 6 . This is 
an attractive property of these peptides that can be used as a scaffold for designing therapeutic peptides, especially 
for improving the stability of linear peptides 711 . In addition, these peptides display cytotoxicity to microbes 12 , 
along with HIV-infected 13 " 15 , cancer 16 and tumor 17 cells. However, their use as therapeutic agents remains limited 
because of their high toxicity to normal cells 8 . Therefore, many efforts have been made to study the mechanism of 
their activity 18 " 26 . 

Cyclotides are amphiphatic peptides with a partitioning between the clusters of hydrophobic and hydrophilic 
amino acid residues on their surface 22 . These clusters are responsible for the membrane binding 19 " 21 and self- 
assembly 22 of the peptides, which are fundamental steps in their activity 18,22 " 24 . Furthermore, it has been proposed 
that cyclotides are able to penetrate the membrane, causing the formation of "barrel-stave" or "toroidal" mem- 
brane pores with an estimated diameter of 4.1-4.7 nm 23,26 . The activity of cyclotides is also involved with the 
depletion of membrane lipids 25 , which is believed to be due to the solubilization of membrane lipids by the 
peptide 26 . Based on those findings, Wang and colleagues 26 have proposed a model to describe the overall mech- 
anism of the activity of kalata Bl (kBl) -a cyclotide in the Mobius subfamily. 

Recently, the molecular dynamics (MD) simulation method, which provides details of molecular motions and 
interactions as a function of time, has been widely used to study the dynamic interactions of biological macro- 
molecules 27,28 . A modern MD simulation technique known as coarse-grained molecular dynamics (CG-MD) 
simulation has been successfully used to study heterogeneous membrane models 29 , which is an updated model of 
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Figure 1 | Progressive membrane association and disruption activity of kBl. The binding of (a) the tower-like and (b) wall-like oligomer of kBl to the 
membrane-water interface, (c) and (d) The top and side view of the ring-like hollow in the kBl cluster at the membrane surface of 120 kBl molecules 
indicating the curvature of the membrane induced by the binding of kB 1 , by which lipids covered by the ring-like hollow were at the same level as the kB 1 
cluster, (e) The side-view snapshot of the simulation of 350 kBl molecules with the membrane indicates that lipids were extracted from the membrane via 
the ring-like hollow at the membrane surface and moved to the space inside the kBl cluster. For all snapshots, the hydrophobic and hydrophilic residues 
of kBl are represented by white and green colors, respectively. DPPC, DUPC and CHOL are colored as lime, cyan and yellow, respectively. Head (NC3) 
and first acyl-chain carbon (CIA and C1B) atoms of DPPC and DUPC are shown as red and blue color, respectively. Molecules of kBl are shown as 
surface models. Lipid molecules are represented as licorice models. To clarify the picture, water molecules are not shown, except (e) that water molecules 
are shown with a violet color to indicate spaces inside the kBl cluster. Each scale bar is equal to 1 nm. 



biological cell membrane 30 . The organization of various peptides 
upon bonding with the membrane model has also been reported 31,32 . 
In addition, the aggregation behavior of cyclic antibacterial peptides, 
which leads to membrane disruption, has also been successfully 
described using CG-MD simulations 33 . Accordingly, we applied 
CG-MD simulations to dynamically investigate the mechanism of 
kBl activity in relation to the membrane disruption. In this study, the 
simulations were performed at various concentrations of kBl 
(Supplementary Table SI online). A simulation of 24 molecules of 
kBl (type 1 simulation) was initially performed to understand the 
progressive oligomerization of kBl. The membrane disruption activ- 
ity of kBl was further observed at higher concentrations of kBl. In 
addition, the effect of the aggregation behavior of kBl on membrane 
disruption activity was investigated in this study. Therefore, the con- 
centration of kBl was increased in two different manners. First, the 
concentration of kBl was increased stepwise from 24 to 48, 72, 96 
and 120 molecules (type 2 simulation). Second, a set of simulations 
with 48, 72, 96 and 120 molecules of kBl were separately performed 
(type 3 simulations). An additional simulation with 350 kBl mole- 
cules was further performed to investigate the membrane lipid 
extraction activity of kBl. Consequently, an alternative view of the 
mechanism of kBl activity that resembles the concept of the inter- 
facial activity model 34 , which is a non-pore membrane disruption 
model 35 , was demonstrated (Fig. 1). 

Results 

Behavior of kBl oligomers. An analytic ultracentrifugation study 
indicated the ability of cyclotides to form oligomers 18 . Another study 
demonstrated that the lateral diffusion of kBl on the membrane 
surface allowed the peptide to form an oligomer containing 16- 
24 kBl molecules that was able to form a ring in the membrane 23 . 
Here, the formation and behavior of kBl oligomers during their 
association with the membrane were examined with a simulation 
of 24 kBl molecules (type 1 simulation). Initially, dimers, tetramers 
and octamers were determined to be the majority of the oligomeric 
forms (Figs. 2a and b). Afterwards, oligomers containing 12, 16 and 
24 kBl molecules were continuously formed (Figs. 2c-e, respectively). 
Once the 24-mer was formed, the kBl molecules did not disassociate 
(Supplementary Fig. SI online). These results indicate that the number 
of kBl molecules in an oligomer was an even number and a factor of 
the tetramer or octamer, which agrees with the results suggested in 



(ref. 23). In addition, the analysis indicated that the oligomers grew 
quickly during 0-0.14 (is, when most of the kBl molecules or 
oligomers were in the water. After that time period (0.15-1.88 (is), 
the oligomers grew slowly when most of the kBl molecules were 
bonded to the membrane. 

For 60 pis during the type 1 simulation, we detected two main 
configurations of kBl oligomers, termed tower-like and wall-like 
clusters (Figs, la and b, respectively). The tower-like oligomers were 
detected during the early period of the simulation, when the kBl 
oligomers in the water bonded to other kBl molecules or oligomers 
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Figure 2 | Clustering behavior of kBl molecules. The size of the kBl 
clusters is plotted as a function of time during the first 2 us of the type 1 
simulation, whereas the number of the clusters in each cluster size is 
represented by the color (see color scale bar). The plot for the full time of 
the simulation is shown in Supplementary Fig. S 1 online. Snapshots for 
each period of time are labeled and covered by the same letter and box color 
as those shown in the plot. Molecules of kBl are shown as surface models 
and are colored by molecular chains. 
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at the membrane surface, resulting in vertically-grown oligomers. 
The tower-like oligomers were then transformed to the wall-like 
oligomers, where most of the kBl molecules were bonded with the 
membrane and longitudinally arranged along the membrane surface. 
When the concentration of kBl was increased (type 2 and 3 simula- 
tions), the wall-like clusters were elongated (Fig. 3). When 120 kBl 
molecules were simulated, the conjugation between the wall-like 
clusters resulted in the formation of a ring-like hollow in the kBl 
cluster ('ring-like' refers to a closed ring made by the connection of 
kBl molecules and 'hollow' refers to space inside the kBl ring, see 
Fig. lc). The hollow was located on the membrane surface because 
kBl interfacially bound to the membrane. The diameter of the ring- 
like hollow detected in the type 2 simulation was approximately 
13 nm (Fig. 3e), whereas the diameter detected in the type 3 simula- 
tions was approximately 5.5 nm (Figs, lc and 3i), which was closer to 
4.1-4.7 nm estimated by Huang et al. (ref. 23). 

Interfacial membrane binding of kBl. The distance between the 
COM of each peptide molecule to that of the membrane was 
measured to determine the position of the peptides relative to the 
membrane. We observed that there were no kBl molecules in any 
simulation that penetrated the membrane (Fig. 4). The average 
distances from the center of mass (COM) of the membrane-bound 
kBl molecules in the type 1 simulation, for example, to that of the 
membrane during 50-60 us of the simulation were in the range of 2.2 
± 0.3(s.d.) and 4.3 ± 0.3(s.d.) nm (with minimum and maximum 
values of 1.4 and 5.2 nm from the raw data, respectively), in which 
the relative distance of the lipids head, glycerol and first acyl- chain 
carbon (CIA and C1B) atoms of the membrane were 1.9 ± O.O(s.d), 
1.4 ± 0.0(s.d.) and 1.1 ± 0.0(s.d.) nm, respectively. These results 
demonstrate that kBl molecules were located at the membrane- 
water interface. 

Membrane deformation by kBl. To investigate the membrane dis- 
ruption activity of kBl, membrane properties changing due to the 
binding of kBl s was observed. Instead of penetrating the membrane, 
the simulations indicated that kBl induced membrane curvature. The 
membrane curvature degree was calculated by measuring the distance 
of every lipid head atom relative to the COM of the membrane. 



Thereby, the positions of all lipid head atoms relative to the COM 
of the membrane were determined (Fig. 4). The results indicated that 
the curvature could be detected even when the membrane was not 
bound with kBl molecules (Fig. 4a). The frequency distribution 
analysis of the distance values revealed two main peaks of the 
distribution graph at 2.0 and —2.0 nm, reflecting the upper and 
lower surfaces of the kBl -free membrane, respectively (Supplemen- 
tary Fig. S2 online). The maximum distances relative to the COM of 
the membrane of this system were 3.2 and —3.5 nm for the positive 
(up-curving) and negative (down-curving) curvature, respectively. 
These results indicate the predominantly negative curvature of the 
kBl-free membrane. However, when the membrane was bound with 
24 molecules of kBl, the maximum distance for the positive and 
negative curvatures, which averaged over 60 us in the type 1 
simulation, increased by 16% (3.7 ± 0.4(s.d.) nm) and 3% ( — 3.6 ± 
0.2(s.d.) nm), respectively. These results indicate that the membrane 
binding of kBl dominantly induced positive curvature. Furthermore, 
the results from the type 1 (Fig. 4b), type 2 (Fig. 4c-f) and type 3 
(Fig. 4g-j) simulations were compared to determine the relationship 
between the concentration of kBl and the degree of the membrane 
curvature. The data in Table 1 indicate that the peak-to-peak 
amplitudes of the membrane curve with 48, 72, 96 and 120 kBl 
molecules in the type 2 simulations increased by approximately 4, 

16, 28 and 49% from that of the type 1 simulation during the same 
period of measured time (12 (is), respectively. In addition, the results 
also indicated that the membrane curvature of the type 2 simulation 
was increased both in the positive and negative direction and spanned 
on a large area of the membrane (Fig. 4 and Supplementary Fig. S2 
online). However, the peak-to-peak amplitudes with 48, 72 and 
96 kBl molecules in the type 3 simulations (Fig. 4g-i) decreased by 

17, 21 and 2% from that of the type 2 simulation at the same 
concentration, respectively. Correspondingly, the data in Table 1 
indicate that the number of kBl molecules that bound with the 
membrane in the type 3 simulations was lower than that observed 
in the type 2 simulation at the same concentration. The results 
demonstrated that the degree of the membrane curvature was 
directly affected by the number of kBl that associate with the 
membrane. Surprisingly, the membrane curvature degree with 
120 kBl molecules in the type 3 simulation (Fig. 4J) was 85% 
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Figure 3 | Different oligomerization behaviors of kBl. (a) The top view snapshot at 12 las of the type 1 simulation indicates short wall-like clusters on the 
membrane surface, (b-e) The wall-like clusters were continuously elongated when the concentration of kBl was increased to 48, 72, 96 and 120 molecules 
in the type 2 simulation, (f-i) The large oligomers were detected on the membrane surface of the type 3 simulations. The representation is the same as in 
Fig. 1, except that the NC3, CIA and C1B atoms are not highlighted. The pink line shown in (e) represents a ring-like hollow in the kBl cluster with a 
diameter of ~ 13 nm, whereas the red line in (i) shows a ring-like hollow in the kBl cluster with a diameter of —5.5 nm. Yellow dashed-line boxes 
represent central simulation cells in a four-periodical image of each simulation system. 
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Figure 4 | Membrane curvature change due to the binding of kBl. The membrane curvatures at various kBl concentrations and simulation times are 
represented by positions of all lipid head (NC3) atoms for both upper and lower leaflets of the membrane relative to the COM of the membrane. The color 
of each lipid head atom is represented by its distance relative to the COM of the membrane (see color scale bar). The kBl molecules are shown as 
transparent surfaces with a violet color on the top of the membrane surface. Top view snapshots of (b-j) are shown in Figs. 3a-i. 



higher than that of the type 1 simulation at 12 ms (Fig. 4b). In 
addition, the value was increased by 25% from that of the type 2 
simulation at the same concentration (Fig. 4f), although the 
number of the membrane-bound kBl molecules in the type 3 
simulations was lower than that observed in the type 2 simulation 
(Table 1). The side-view snapshot of the ring- like hollow in the kBl 
cluster in the type 3 simulations indicates a relationship between the 
ring-like hollow formation and the degree of the membrane curvature 
because it showed that the membrane in the area that surrounded by 
the ring-like hollow (which was the highest membrane curve among 
all systems) was at the same height level as the kBl molecules that 
forming the ring-like hollow (Fig. Id). Quantitatively, the relative 
position between the membrane and kBl cluster was represented 
by their time series hydrophobicity ratio (Fig. 5), which was 
measured perpendicularly to the membrane surface (see Supplemen- 
tary Fig. S3 online for details of the measurement). The hydropho- 
bicity ratio can also represent thickness change of the membrane. 
Initially, the thickness of the membrane was ca. 5 nm (see Fig. 5a), 
where the lower and upper surface of the membrane were located at 
ca. 5 and 10 nm along the Z axis of the simulation box, respectively. 
The hydrophobicity ratio at the membrane core was in range of 0.8-1 
while those of upper and lower hydrophilic layers were in range of 
0.6-0.1. The membrane thickness was then rapidly increased to 
approximately 16 nm. This thickness change was a result of the 



membrane curvature, as indicated by the lower hydrophobicity 
ratio at the center of the thickness-increasing area of the mem- 
brane (see dashed-line box) compared with that at the core of the 
membrane at 0 ms. Interestingly, 85% of the total increase of the 
thickness occurred in the upper layer, whereas only 15% was 
detected in the lower layer of the membrane, indicating a dominant 
positive curvature over the negative curvature. For the cluster of kBl 
molecules of the same system as the membrane (Fig. 5b), the 
hydrophobicity ratio was high at the bottom region of the cluster 
(see solid-line box). The kBl molecules in this region bound with 
the membrane. However, a low hydrophobicity ratio was also 
detected at the lowest area of this region because some kBl 
molecules in this area could not bind to the membrane due to the 
membrane curvature (Fig. Id for snapshot). A low hydrophobicity 
ratio was also detected at the top of the kBl cluster, which was the 
region that kBl molecules interacted with the water. In contrast, the 
hydrophobicity ratio at the core of the kBl cluster was in range of 
0.5-0.6. In the beginning of the simulation (the curvature was not 
present in the membrane) the kBl cluster was on top of the 
membrane and its hydrophobicity-plot area was slightly overlapped 
with that of the membrane (see Fig. 5a and b during 0-0.1 u.s at ca. 8- 
9 nm of the Z axis). This corresponds with the interfacial membrane 
binding of kBl. Interestingly, when the membrane curvature was 
induced by kBl, the thickness- increasing area of the membrane (ca. 



Table 1 | Summary of membrane curvature changes of the three simulation types related to kBl concentration. 'Max' refers to the maximum 
points of the positive membrane curvature where 'Min' refers to the minimum points of the negative membrane curvature. 'Peak- to- Peak 
amplitude' is the distance between Max and Min points. Values of 'The number of membrane-bound kBl ', 'Max' and 'Min' are a time- 
averaged values, and the error estimation is shown as standard deviation 



Simulation type Time (us) 


Total number of kBl The number of membrane-bound 

(molecule) kBl (molecule) Max (nm) 


Min (nm) 


Peak-to-Peak 
amplitude (nm) 


Type 1 


0- 


-12 


24 


9 ± 1 


3.6 ± 0.3 


-3.7 ± 0.2 


7.3 


12- 


-24 


24 


8 ± 1 


3.9 ± 0.4 


-3.6 ± 0.2 


7.5 




24- 


-36 


24 


1 1 ± 1 


4.0 ± 0.4 


-3.7 ±0.2 


7.7 




36- 


-48 


24 


8 ± 1 


3.7 ± 0.3 


-3.5 ± 0.2 


7.2 




48- 


-60 


24 


8 ± 1 


3.4 ± 0.3 


-3.4 ± 0.2 


6.8 


Type 2 


12- 


-24 


(24 x 2) 


21 ± 1 


3.9 ± 0.3 


-3.9 ± 0.2 


7.8 




24- 


-36 


(24 x 3) 


25 ± 3 


4.8 ± 0.3 


-4.1 ± 0.3 


8.9 




36- 


-48 


(24 x 4) 


28 ± 2 


4.9 ± 0.2 


-4.3 ± 0.2 


9.2 




48- 


-60 


(24 x 5) 


32 ± 2 


5.2 ± 0.3 


-4.9 ± 0.2 


10.1 


Type 3 


0- 


-12 


48 


10 ± 1 


3.2 ± 0.2 


-3.3 ± 0.2 


6.5 


0- 


-12 


72 


16 ± 1 


3.4 ± 0.2 


-3.6 ± 0.2 


7 




0- 


-12 


96 


19 ± 2 


4.6 ± 0.2 


-4.4 ± 0.2 


9 




0- 


-12 


120 


23 ± 2 


6.8 ± 0.2 


-5.8 ± 0.2 


12.6 
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Figure 5 | The membrane was lifted to the kBl cluster. The 

hydrophobicity ratio of (a) the membrane and (b) a cluster of 120 kBl 
molecules located on the membrane surface were measured along the Z 
axis of the system (see Supplementary Fig. S3 online). A dashed-line box 
covers the hydrophobicity ratio at the thickness-increasing area of the 
membrane. A solid-line box covers hydrophobicity ratio at the bottom 
region of the kB 1 cluster. The Y axis of the heat map represents overall 
thickness (Z axis) of the simulation box. 

10-16 nm in the Z axis) significantly overlapped with the area of the 
kBl cluster, whereas the position of the kBl cluster was slightly higher 
during 0-3 (is. These results clearly indicate that kBl molecules did 
not deeply penetrate the membrane. Instead, the membrane was 
curved up onto the kBl cluster. 

Membrane lipid extraction by kBl. An additional simulation was 
further performed with 350 molecules of kBl according to the kBl 
molecule per lipid ratio of 0.35 that caused 50% lipid vesicle leakage, 
as calculated from the graph showing the relationship between the 
ratio of kBl with membrane lipids and bioactivity of the peptide 23 . 
Membrane behavior similar to that previously described was also 
detected in this system. Additionally, the membrane lipid 
extraction was further analyzed by comparing the position of lipid 
residue 296 (as an example) with the COM of the membrane and that 
of the membrane surface. Positions of sample lipid, the whole 
membrane, and the membrane surface were obtained by measured 
distance from their COM to the reference point (see Supplementary 
Fig. S4a online). The results showed that the distances from the COM 
of the whole membrane and that of the membrane surface to the 
reference point were continuously increased as a result of the positive 
curvature. Simultaneously, the distance of lipid residue 296 relative 
to the reference point was also increased (Fig. 6a and b). The distance 
increase of the COM of lipid residue 296 was steady during 1-12 [is. 
Therefore, to observe the increasing rate of the distance, the distance 
values of lipid residue 296, the whole membrane, and the membrane 
surface were fitted using a linear equation during this time period. 
The results indicated that the slope of the linear fitting curve for lipid 
residue 296 (0.27) was approximately 2 and 3-fold higher than that of 
the whole membrane (0.14) and that of the membrane surface (0.10), 
respectively (Fig. 6a). Similar results were also observed in other 
sample lipids (Fig. 6b and Supplementary Fig. S4b online), 
implying that lipid residue 296 and other lipids were extracted 
from the membrane. 

Discussion 

Although the mechanism of kBl activity has been extensively studied 
and several models have been proposed to depict the overview of the 
mechanism, the intensive dynamic and molecular details of the 
mechanism that could provide a deeper insight have not yet been 
reported. Here, CG-MD simulations were conducted to demonstrate 
progressive membrane association, oligomerization and membrane 
disruption activity of kBl. Based on the results from the MD simula- 
tions, we propose an alternative view of the mechanism of kBl activ- 
ity. Fig. 1 shows the overall membrane disruption progression of the 
newly proposed model. The membrane disruption activity of kBl 



molecules is initiated by their binding to the membrane. We 
observed that some kBl molecules directly bound to the membrane, 
and some other molecules formed oligomers in the water. In addi- 
tion, the results suggest that the tetramer was a building block of the 
oligomer. The kBl oligomers in the water then bound to the mem- 
brane surface, forming tower-like clusters (Fig. la). During this per- 
iod (0-0.14 us), the oligomers quickly grew compared with the later 
period (0.15-1.88 us) when most of kBl molecules bound to the 
membrane. This result was likely because kBl in the water could 
move and rotate in all directions to interact with other kBl mole- 
cules. In contrast, the membrane-bound kBl molecules could only 
move on the plane of the membrane surface, thus reducing their 
probability of interacting with others. However, the hydrophobic 
residues of kBl molecules that formed the towers were mostly 
exposed to the water (Fig. la). The tower-like clusters were then 
transformed to the wall-like cluster to bury the hydrophobic residues 
in the membrane (Fig. lb). When the concentration of kBl was 
increased, the wall-like cluster was then elongated by connecting 
with other wall clusters and eventually forming the ring-like hollow 
in the kBl cluster on the membrane surface (Fig. 3). However, we 
observed that a different manner of kBl concentration increasing 
could affect the ring-like hollow formation. In the type 2 simulation, 
the kBl molecules that were stepwise increased formed small oligo- 
mers) in the water before binding to the membrane and then lat- 
erally moved to interact with other oligomers (Figs. 3a-e). This 
allowed for the scattering of kBl molecules on the membrane surface 
and resulted in the formation of the large ring-like hollow in the kBl 
cluster (Fig. 3e). In contrast, kBl molecules in the type 3 simulations 
reached the membrane as large oligomers. Therefore, the smaller 
ring could be formed more easily in this simulation type because 
the kBl molecules were less scattered in this system, and thus, a 
long-distance lateral movement of kBl was not required (Figs. 3f-i). 

Wang and colleagues proposed that kBl molecules are able to 
permeate the membrane to form pores in the membrane 26 . 
Nevertheless, the dynamics and molecular details of the kBl -mem- 
brane interaction obtained in this study demonstrated that kBl 
molecules did not penetrate the membrane and were trapped at 
the membrane surface. We observed that kBl molecules were located 
at the membrane interfacial zone, which is the region at the mem- 
brane-water interface and has a thickness of 1.0-1.5 nm covering the 
water area, where the membrane is polar, with small portions of 
hydrophobic atoms 34 . Therefore, the insertion of kBl in this zone 
could be described by the partitioning of hydrophobic and hydro- 
philic patches on the peptide surface. We performed a short simu- 
lation of one kBl molecule with the membrane to verify this 
assumption and found that with the membrane-bound state, 
approximately 55.84% and 44.16% of the total interaction energy 
of the kBl monomer accounted for the kBl -membrane and kBl- 
water interaction energy, respectively (Supplementary Fig. S5 
online). These results indicate that although the hydrophobic patch 
was inserted into the membrane, nearly half of its structure, which 
was the hydrophilic patch, preferentially interacted with the water. In 
addition, our study used the same simulation force field (MARTINI 
CG force field 36 3S ) as that used in (ref. 39), which demonstrated that 
two of the three forms of cardiotoxin A3, the Cobra Cytotoxin, which 
has a configuration similar to that of cyclotides, lack membrane 
penetration ability while another form is able to permeate the mem- 
brane 39 . Therefore, the lacking of membrane penetration ability of 
kBl should not be a result of the use of MARTINI CG force field. 

Instead of penetrating the membrane, kBl molecules began dis- 
rupting the membrane by inducing curvature in the membrane while 
they were still on the membrane surface. In fact, the membrane 
curvature could be detected even when the membrane was not bound 
with kBl molecules (Fig. 4a). Negative curvature was observed at the 
CHOL-enriched domain of the kBl-free membrane, which is con- 
sistent with the fact that the existence of small lipids, such as choles- 
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Figure 6 | Membrane lipid extraction by kBl. (a) The distance of the COM of a sample lipid (lipid residue 296) and the COM of the whole membrane, as 
well as that of the membrane surface to the reference point, are plotted as a function of time. Green, red and blue lines show the distance from the COM of 
the whole membrane, the membrane surface and lipid residue 296 to the reference point, respectively. Pink, yellow and cyan lines show the linear fitted 
curve of the distance values of the whole membrane, the membrane surface and lipid residue 296, respectively, (b) NC3 atoms of DPPC and DUPC in 
upper and lower layers of membrane are shown as green balls. Lipid residue 296, 496, 585, 629, 748 and 927 are shown as a CPK model with red, 
violet, black, blue, yellow and orange colors, respectively. Size of each grid space is 1 X 1 nm. (c) The volume matrix slides in the X, Y and Z planes of the 
system show the cross-section surface inside the kBl cluster. The areas occupied by kBl molecules are shown as small dots on the volumetric slides, 
whereas the spaces inside kBl cluster are indicated by the arrows. Lipid residue 296 is shown as a licorice model with a red color. Other lipids in the 
membrane are also shown as licorice models with the same color as shown in Fig. 3. 



terol, at some local areas on the membrane could result in negative 
membrane curvature at that area 40 . In contrast, when the membrane 
was bound with kBl, positive curvature was detected in the majority 
areas. The curvature degree (height of the curve) was directly 
dependent on the number of kBl molecules that interfacially bound 
with the membrane (Table 1). These results are consistent with 
reports (ref. 40, 41) that demonstrated that the asymmetric mem- 
brane binding of amphipathic peptides with high local density can 
result in membrane curvature. Furthermore, the insertion of peptides 
into approximately 16% of the total thickness of the membrane 
hydrophobic region resulted in a positive curvature, whereas a nega- 
tive curvature occurred when the insertion depth was 42% of the 
thickness 42 . We then measured the distance of the COM of kBl to 
that of the membrane to explore insertion dept of kBl in membrane 
and observed the distance of 2.2 ± 0.3(s.d.) which was the area 
occupied by the lipid head atoms (Supplementary Fig. S6 online). 
We also calculated the distance of Trpl9, which is the largest 
residue of kBl that inserted into the membrane, to the COM of the 



membrane in order to measure the deepest distance that a kBl mole- 
cule can insert into membrane. The results indicate that the average 
insertion depth of Trpl9 of single kBl molecule was 1.7 ± 0.4(s.d.) 
nm relative to the COM of the membrane (Supplementary Fig. S6 
online). Thus, Trpl9 was inserted into approximately 7% of the 
hydrophobic region, as calculated from the lower bound of the aver- 
age value of 1.3 nm when GL atoms were considered as the surface of 
the membrane hydrophobic region. These results indicate that the 
partial insertion of kBl into the membrane interfacial zone was 
responsible for the membrane curvature, especially the positive 
curvature. 

However, the membrane curvature was more dominantly related 
to the formation of the ring-like hollow in the kBl cluster of kBl at 
the membrane surface. In the type 3 simulations at 120 molecules of 
kBl, we observed a ring-like hollow with a diameter nearly identical 
to that reported with the experiment 23 (Fig. lc and 3i). We also 
observed that the membrane in the area that surrounded by the 
ring-like hollow (which was the highest curvature area among all 
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simulation systems) was curved to the same level as the kBl that 
formed the ring-like hollow (Fig. Id). This relative position between 
kBl molecules and the membrane, i.e., position of kBl molecules 
were lower than the upper surface of the membrane (see results and 
Fig. 5), is similar to the results reported in (ref. 26). However, the 
movement of the membrane was limited in (ref. 26) because it was 
held on the support surface according to the neutron reflectometry 
method used in that study, and it was proposed that the cyclotides 
completely inserted and penetrated into the membrane 26 . In contrast, 
both the kBl and membrane in our study were freely dynamic, and 
we molecularly demonstrated that kBl did not penetrate the mem- 
brane; it was the membrane that curved up onto the kBl cluster. 

Burman and colleagues demonstrated that cy02 in the Bracelet 
subfamily of cyclotides was able to extract lipids from the mem- 
brane 25 . Our results consistently demonstrated that kBl, even if it 
is a member of different subfamily, also displayed membrane lipid 
extraction activity. Our results indicated that kBl molecules formed a 
large cluster on the membrane surface. Moreover, an analysis of kBl 
occupancy demonstrated that there were spherical spaces occupied 
by water molecules inside the kBl cluster and were connected to the 
ring-like hollows at the membrane surface (Fig. 6c). The overlaying 
of the volumetric slide of the kBl occupancy to the membrane sur- 
face indicated that the lipids were extracted from the membrane 
through the kBl ring-like hollow at the membrane surface and then 
moved into spaces inside the kBl cluster (Fig. 6b and c). Therefore, it 
can be said that the ring-like hollow in the kBl cluster was a starting 
area of membrane disruption by kBl. The movement of lipids from 
the membrane to the kBl cluster was possibly involved with hydro- 
phobicity distribution in the kBl cluster. The ratios between hydro- 
phobic and hydrophilic atoms at the core of the kBl cluster were 0.53 
± 0.03(s.d.) and 0.47 ± 0.03(s.d.), respectively (Fig. 5b). This might 
have resulted in stepwise movement of lipids in the kBl cluster by 
which acyl-chains of the lipids were attracted by hydrophobic resi- 
dues while lipid head groups were attracted by hydrophilic residues. 
However, the hydrophobic mismatch between the lipids and kBl 
molecules in the cluster could have occurred because the hydropho- 
bicity ratio in the oligomers was not substantially different. Thus, 
there could have been continuous movement of lipids due to the 
chain of favorable and unfavorable hydrophobic interactions 
(Fig. 6b). 

In summary, the mechanism obtained in this study reflects an 
interfacial activity model that used to describe the activity of anti- 
microbial peptides 34,35 . The interfacial activity model describes the 
ability of a molecule to bind to the membrane, partition into the 
membrane-water interface, and alter the packing and organization 
of the lipids 34 . In addition, this model was characterized as a non- 
pore model for the membrane disruption activity of antimicrobial 
peptides 35 . Thus, our study provides an alternative view of the mech- 
anism of kBl activity. 

Methods 

All of the simulations were conducted with the GROMACS software package 43 ' 44 
using the MARTINI CG force field 36-38 . The CG model is a molecular representation 
technique that reduces the number of atoms in the system by mapping a set of atoms 
into pseudo-atoms 37 and is thus suitable for long simulations and complex processes 
of pep tide-membrane interactions 45 . The simulation integration was performed using 
a leap frog algorithm with a time step of 20 femtoseconds. The temperature was 
independently coupled for each molecule in the system to 3 10 K using the Berendsen 
thermostat 46 . The pressure was coupled to 1 atm using the semi-isotropic scheme of 
the Berendsen barostat 46 , where the x, y plane and the z direction were coupled 
separately to obtain a tensionless membrane. The periodic boundary condition with 
standard non-bonded interaction criteria was applied to all systems. Visual Molecular 
Dynamics (VMD) software version 1.8. 6 47 was used for system visualization and 
analysis. Because the simulations were performed based on a CG model, the time 
shown in the results was scaled by a factor of 4 (see ref. 29 for details). 

Preparation of the KBl CG model. Atomic coordinates of KBl were obtained from 
the Protein Data Bank (PDB code 1NB1) 3 . The file "atomlCg_v2.1.awk" was 
downloaded from the MARTINI website and used to convert the atomic coordinates 
to a CG coordinate file. Most of the amino acid residues were represented by two 



pseudo-atoms, one for the main chain and one for the side chain. Tryptophan was 
represented by five pseudo-atoms, one for the main chain and four for the aromatic 
side chain. Arginine was represented by three pseudo-atoms, one for the main chain 
and two for the longer side chain, whereas glycine was represented by one pseudo- 
atom showing only the main chain. The amino acid sequence and the secondary 
structure information were used as the input of "seq2itp.pl" to generate the 
GROMACS topology file. The sequence was obtained from the Protein Data Bank, 
and the secondary structure information was prepared using the GROMACS utility 
program do_dssp. Two additional options of "cysteines" and "elasticin" from the 
seq2itp.pl script were used to automatically define the disulfide bond and local elastic 
networks for kBl. Because the N and C termini of kBl are closed by a peptide bond, 
the generated topology file was further modified to reflect this. Based on the sequence 
of 1NB1, in which the N terminus is cysteine and the C terminus is valine, all bond 
parameters of Val2 1 connecting with Cys 22 were used to describe the peptide bond of 
Val29 with Cysl. In this work, all simulations were performed at a pH equal to 7, and 
thus, the net charge of kBl was 0. 

Preparation of the membrane model. The heterogeneous membrane model was 
constructed based on lipids ratio obtained from experimental study 48 in order to 
mimic the real cell membrane. The number of lipids model was calculated based on 
the percent of lipid composition of uninfected H9 cell membrane. For 100% of total 
phospholipids, approximately 9.8% were sphingolipids, and 90.2% were other 
phospholipids 48 . The cholesterol composition was expressed as the cholesterol- to - 
phospholipid (c/p) ratio of 0.476. Using CG- MARTINI force field 36 " 38 , we applied the 
lipid composition regarding (ref. 29), including diundecanoyl-phosphatidylcholine 
(DUPC), dip almitoyl- phosphatidylcholine (DPPC), represent saturated 
(sphingolipid) and unsaturated lipids, respectively, and cholesterol (CHOL). In our 
system, there are 1,000 lipid molecules with the DPPC : DUPC : CHOL ratio of 
0.07 : 0.62 : 0.3 1 . in total. The lipid system was hydrated by six water beads (equal to 24 
real water molecules) per lipid molecule 49 . Structures and GROMACS topology files 
of these three lipids were obtained from the martini website (http://md.chem.rug.nl/ 
cgmartini/). We started construct the membrane model by uploading single molecule 
of DUPC, DPPC and CHOL to the web service of our in-house software, Automated 
Membrane Generator (https://10.4.52.26/amg/). The program automatically 
generated lipid bilayer of DPPC and CHOL with the ratio of 1 : 1 based on the 
theoretical shape complementary described in (ref. 50). The domain was then 
surrounded by DUPC and the remaining CHOL. The output heterogeneous lipid 
bilayer was downloaded from the web service. The lipids ratio in the bilayer was 
modified to obtain the desired membrane. After energy minimization, a simulation of 
100 ns was performed for system equilibration at a temperature of 310 K and 
pressure of 1 atm. The production run of the 4 \is simulations was conducted under 
the same conditions. The average DPPC : DUPC : CHOL ratio in the domain area of 
the constructed membrane model was about 0.53 ± 0.20(s.d.) : 0.07 ± 0.17(s.d.) : 0.39 
± 0.23(s.d.), which is nearly identical to that reported in (ref. 29). The thickness 
difference between the domain and non-domain area of the membrane was in the 
range of 0.4-0.7 nanometers, which was corresponded to that of previously observed 
in (refs. 29,5 1 ,52). To validate our membrane model, we increased the temperature up 
to 400 K to disorganize the lipid domain formation according (ref. 29). The 
temperature was then cooled down to 310 K. After energy minimization and 
equilibration, the production run of this system was performed for 20 us. At this time 
the domain was reformed (Supplementary Fig. S7 online). 

Simulations to investigate oligomerization and membrane disruption activity of 
kBl. Three simulation types relating to kBl concentration were performed to 
investigate the progressive membrane binding, oligomerization and membrane 
disruption activity of kBl. Because the charge distribution on the surface of kBl is 
asymmetrical, the molecules of the peptides were added into the system with random 
geometry orientation. The molecules of kBl were randomly added into the system 
using the Genbox program (see http://manual.gromacs.org/current/online/genbox. 
html) close to the upper surface of the membrane. All simulation types were started 
with the same membrane, the membrane at 4 us, which contained 1,000 lipids 
molecules and ca. 24,000 water molecules (equal to 96,000 real water molecules). For 
the type 2 simulation, the kBl molecules were continuously added into the system 
every 12 us by replacing water molecules. Then, the water molecules were added for 
each kBl concentration to meet the desired amount of water molecules. Therefore, 
thickness and size of the systems increased when the number of kBl molecules was 
increased (see Supplementary Table SI online). In addition, we placed the membrane 
at about 10-20% lower than center of the simulation box in the Z axis and then placed 
the peptides on the top of the membrane. 

Size analysis of the kB 1 oligomer. The clustering of kB 1 in the oligomer was 
analyzed using the GROMACS utility tool-g_clustsize (http://manual.gromacs.org/ 
online/g_clustsize.html). The "-mol" option of the program was used for calculating 
the clustering of the molecules. The cut-off for the clustering (see "-cut" option) was 
2.0, which we determined was the average distance between two kBl molecules when 
bound together. The default settings for all other options of the program were used. 

Analysis of hydrophobicity ratio. The system was divided into small slabs along the 
Z axis. The numbers of hydrophobic and hydrophilic atoms in each slab were counted 
for the membrane and kBl clusters. The hydrophobic ratio was calculated by dividing 
the number of hydrophobic atoms by the total number atoms in each slab. The 
hydrophobic ratios in all slabs were combined together based on their position along 
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the Z axis. The process was repeated for every time step of the simulation (200 ps) and 
performed separately for kB 1 and the membrane. This analysis was performed using a 
custom tcl-script implemented by using the VMD 47 library. 

Analysis of the occupancy volume of kB 1 . The occupancy volume was calculated 
using the VolMap Plugin of the VMD 47 software (http://www.ks.uiuc.edu/Research/ 
vmd/plugins/volmapgui/). The analysis was performed for all residues of all kBl 
molecules. The volumetric data of the "map types" option was set to "density", with 
weighting on the occupancy. The values were set to 0.5 and 1.5 A for the "resolution" 
and "atom size" options, respectively. The occupancy volume was averaged for all 
frames of the simulation. The volume was visualized using VMD "VolumeSlice" 
drawing method in the X, Y and Z axes of the system. 

Analysis of the exposure of kBl residues to water. The exposure of the kBl residues 
was represented by the solvent accessible surface area (SASA), which was calculated 
for the hydrophobic, hydrophilic and all other residues of kB 1 . The SASAs were 
analyzed using the "g_sas" program of GROMACS. The solvent probe radius was set 
to 0.56 nm following the method used in (ref. 53). 
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